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PACS 85 . 25 . Dq - Superconducting quantum interference devices (SQUIDs) 
PACS 41 . 60 . -m - Radiation by moving charges 
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Abstract. - We consider the superconducting phase in a moving superconductor and show that 
it depends on the displacement flux. Generalized constitutive relations between the phase of a 
superconducting interference device (SQUID) and the position of the oscillating part of its loop 
are then established. In particular, we show that the Josephson current and voltage depend 
on both the oscillator position and velocity. The two proposed relativistic corrections to the 
Josephson relations come from the macroscopic displacement of a quantum condensate according 
to the (non-)Galilean covariance of the Schrodinger equation, and the kinematic displacement of 
the quasi-classical interfering path. In particular, we propose an alternative demonstration for the 
London rotating superconductor effect (also known as the London moment) using the covariance 
properties of the Schrodinger equation. As an illustration, we show how these electromotive effects 
may induce self-sustained oscillations of a mechanical system. 



Introduction. — The interest to nanomechanical sys- 
tems dramatically increased recently. Last year, the su- 
perposition of quantum states of a mechanical resonator 
has been demonstrated for the first time [Ullj, realizing 
an important step towards the generation of mechanically 
quantum dressed state at the mesoscopic level. These 
dressed states open wide possibilities for using mechanical 
resonators for quantum-limited detection and for quan- 
tum information. There are currently several routes being 
explored towards these applications, based on coupling of 
mechanical resonators to electrical, optical, or magnetic 
systems. One of the routes, theoretically proposed is to 
use a mechanical resonator embedded into a superconduct- 
ing quantum interference device (SQUID), see e.g. [3] and 
references therein. The first experiment demonstrated the 
possibility of the detection of the resonance frequency and 
the quality factor of the resonator by measurements of 
the voltage generated across the SQUID [3]. Magnetic 
flux through the SQUID and the bias current served as 
two control parameters. Later experiments [S] investigated 
back-action of the SQUID exerted on the mechanical res- 
onator and found qualitative agreement with the results 
of the theoretical modeling. 



So far, the description of the experimental setup ignored 
possible electromotive effects of quantum nature. For in- 
stance, it is well known that a cold quantum gas exhibits 
vortex states under rotation, in a intuitive analogy with 
the Abrikosov lattice [S]- Nevertheless, the analogy is un- 
correct, because the Abrikosov lattice is of electromagnetic 
nature, whereas the circulation vortex lattice is of mechan- 
ical origin. A situation when both these effects may com- 
pete is precisely the situation when a quantum condensate 
made of charged particles is mechanically displaced. Then, 
a superconducting condensate may magnify electromotive 
effects when put under displacement. This is illustrated by 
the striking Meissner or London momentum effects [71[5]. 
In short, the Meissner effect corresponds to the genera- 
tion of a displacement current which screens an applied 
magnetic field, whereas the London momentum effect cor- 
responds to the generation of a macroscopic magnetic field 
which screens the displacement current generated by the 
mechanical rotation of a superconductor. This surprising 
effect can only exist when electromagnetism and mechani- 
cal displacement compete together, and can be seen as the 
destruction of the mechanically induced vortex lattice by 
the generation of an electromagnetically induced lattice 
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Another explanation of this effect lies in the well-known 
London theory of electrodynamics of superconductors, 
which corresponds to the addition of an inertial term (pro- 
portional to the vector potential A oc j) into the other- 
wise viscous expressions for circuit electromagnetism (i.e. 
Ohm's theory). Because of this inertial correction, the 
London theory is unable to take into account a change 
of the inertial frame, in the sense that there is no explicit 
need to specify in which inertial frame the superconductor 
is supposed to be. Nevertheless, when a normal metal is 
attached to a superconducting one, the normal electronic 
flow must be recovered at the interface. Then, London 
proposed to correct his theory by imposing a supercon- 
ducting current to lag behind the lattice one, creating a 
magnetic moment by virtue of the Ampere law. To test 
the validity of this retarded contribution, London designed 
the very simple experiment of the rotating superconduc- 
tor, which predicts the generation of a macroscopic mag- 
netic field induced by the rotation of a superconducting 
sphere [7]. The London's prediction was soon after veri- 
fied for both bulk [TU] and proximity effect [TI] systems, 
and latter for high-temperature superconductors [12j in- 
cluding heavy- fermion compounds |13j . 

More recently, the London expression for the inertial 
current was considered in the framework of an effective 
theory of elastic superconductors, with prediction of a 
precise acoustic sensor using Josephson systems, elasto- 
magnetic coupling between the motion of the supercon- 
ductor and the internal magnetic moment it produces, in 
addition to some interesting effects in type-II supercon- 
ductors, see e.g. [M] and references therein. It also con- 
tinue to attract some fundamental interests, being at the 
cornerstone of mechanics and electromagnetism |15| . 

In this paper, we consider some possible electromotive 
effects of quantum mechanical origin in a moving super- 
conductor. We will concentrate on quasi-classical equa- 
tions of motion, and first consider the problem of the 
renormalization of the Josephson relations in a moving 
superconducting system in general terms. Then, using 
simple models and arguments, we will derive constitutive 
relations linking the oscillator motion to the current and 
voltage, and propose an alternative demonstration for the 
London momentum effect based on the covariance of the 
Schrodinger equation. Subsequently, we apply the argu- 
ments to the setup of a suspended SQUID. We believe, 
however, that much of this work can be adapted to other 
detection schemes, especially those using superconducting 
cavities [16] or cold quantum gases [17] . As a simple illus- 
tration of the theory, we will describe the regime when a 
static current can induce self-sustained oscillations of the 
nanoscale bar embedded into the superconducting loop. 

Model of a suspended SQUID. Our starting 
point is the setup shown in FiglTJa which consists of a 
superconducting loop with two Josephson junctions, cou- 
pled to a mechanical harmonic oscillator. The SQUID is 
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Fig. 1: The scheme of a nanomechanical oscillator embedded 
into a SQUID circuit. The assumption of this work is to con- 
sider the scheme a) as being equivalent to the scheme b) (for 
representation simplicity, we drop the self-inductance L, capac- 
itance C and resistance R on the 6) subfigure). In particular, 
this assumption allows us to disregard the elastic properties of 
the bar, and to use relativistic arguments when the bar and 
circuit are two different sub-systems of the entire circuit. 



characterized by two degrees of freedom, which we take to 
be the sum and the difference of the phases of the junc- 
tions, respectively, ip±. The equations of motion can be 
easily written in the quasi-classical approximation. 



i — sin(/)+ cosi^_ + io^ h>+ + [RC /lOc) 
j ~ sin iy9_ cos if+ + ojjT"'-!^- + {RC/u!c) ifi- 



(1) 



with ujc = 27ri?/o/$o the characteristic frequency of the 
SQUID, Iq the critical current of each Josephson junction, 
$0 ~ Trh/e the superconducting flux quantum, R and C 
the resistance and capacitance of the shunted Josephson 
junctions, i = I/2Iq and j ~ J/2/o, where / and J are 
the bias and self-circulating current, respectively [H]. The 
overdot refers to global time derivative. 

When one arm of the superconducting loop is oscillating 
(see FigUJa), Eqs.(|T]) are supplemented by the equation 
of motion for the resonator with the mass m, the length 
£, the quality factor Q, and the resonance frequency wq 
according to 



^0 . , 2 
-X + u,. 



5 (« + j) 



(2) 



where g ss £BIo/m is a geometric acceleration, such that 
g (i + j) represents the Lorcntz force when B is the ex- 
ternal magnetic field. Usually, the SQUID and oscillator 
dynamics are coupled through constitutive relations for 
the phase difference and the phase sum [5H51[T5] 



K^x — f3j and tp-f 



$0 
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-V 
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where /? = AttLIq/^o is the self-inductance of the loop, 
K_ w ttB£/^q is the flux-to-phase geometric coupling. 
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(/?e = TT^/^Oi ^ the external flux across the loop. The 
phase sum-voltage relation has been used to obtain the 
system ([T]) in a quasi-classical way. 

Let us shorty address an illustrating situation in the 
static case. Taking (3 = for simplicity, the phase dif- 
ference If- degree of freedom is frozen by the mechanical 
displacement, and the dimensionless static potential of the 
above situation is 

u {ip+ , x) = 1 — cos ip^ cos {ipe + Kx) + iip+ + oj'^x'^ /2 (4) 

where = uiloq/Ej is the rescaled frequency, and Ej 
is the Josephson energy. The minima of the potential 
correspond to the possibility for the global system to be in 
a fluxon state corresponding to a given amplitude of the 
mechanical bar. Then, it becomes possible to change the 
amplitude of the mechanical deformation and the fluxon 
state by the application of the current i or the magnetic 
flux (fie, and vice-versa. In other words, when a mechanical 
arm of the SQUID loop is able to move, the position x 
degree of freedom is able to affect the magnetic degree of 
freedom of the entire SQUID. Reciprocally, a change in 
the fluxon state of the SQUID can change the equilibrium 
amplitude of the mechanical arm. The amplitude degree 
of freedom is thus formally identical to an electromagnetic 
degree of freedom. 

Alternatively, Eq. ([2]), when seen as the equation of mo- 
tion for the X degree of freedom of the circuit, can be inter- 
preted as the harmonic generation of electromagnetic field. 
Moreover, it is well known that electromagnetism is such 
defined that both position x and velocity x play the role of 
sources [19j . The following discussion aims at introducing 
velocity dependencies into the constitutive relations ([Sj 
and the equations of motion ([I} in order to overcome the 
conceptual difficulties of dealing with mechanical displace- 
ment and electromagnetism in superconducting circuitry. 
Our subsequent reasoning will transform Eqs.((T]) and Q 
into Eqs.Jlil) and (fT5| . 

In order to find the electromotive contributions, we as- 
sume that the elastic scheme described in FiglT]a can be 
made formally identical to the one depicted on FiglTJfe. In 
this way, we transform the elastic bar into a rigid one, 
which can move harmonically on top of the remaining cir- 
cuit. This assumption transforms the previous compli- 
cated situation into a simpler one, quite familiar when 
discussing electromotive effects in term of classical elec- 
tromagnetism [191. In particular, this model allows us to 
make the following assumptions, 

(i) The oscillator displacement function x (t) is a function 
of the time only. 

(ii) The two systems of the oscillating bar and the re- 
maining circuit are relativistically separable. 

These two assumptions are crucial for the following dis- 
cussion. Clearly, without the first assumption, the dis- 
placement function would be a true field dependent on 



both position and time, which requires a complete descrip- 
tion in term of quasi-classical elasto-electro-magnetism. 
This rather complicated situation is out of the scope of 
the present study. Thus, assumption [i] allows us to dis- 
regard some possible elastic ambiguities of the problem. 
The second assumption is just a convenient one in order 
to verify that the forthcoming electromotive terms can be 
justified on the basis of Galilean relativity, eventually com- 
plemented by non-inertial effects. In particular, assump- 
tion (ii) allows us to attach an inertial frame of reference 
to the circuit at rest, and a non-inertial frame to the os- 
cillating bar. 

In order to complete the circuit, the contact points Ki 
and K2 on FiglT]6 are assumed to be frictionless, and elec- 
trical contacts are taken to be ideal. 

Electromotive contributions. Equations ([T|) and 
^ do not take into account the specific properties of a 
moving superconductor, which are expressed as relativis- 
tic contributions to the phase-flux and phase-voltage con- 
stitutive relations ([3]). In order to obtain these contri- 
butions, we constraint the phase of the superconducting, 
macroscopic wave function to be continuous all around the 
circuit loop. When the loop is pierced by a magnetic flux 
^ = § A.dl, the phase continuity imposes to consider the 
so-called gauge covariant phase 7 = <po — 27r$/$0; where 
(/9o is the initial condition phase |18) . 

In the situation depicted in FiglT]6, this gauge covari- 
ant phase fails to describe the passage from the circuit 
at rest to the moving bar. Indeed, 7 is the gauge covari- 
ant phase for circuits at rest. From general properties 
of quantum mechanics, the displacement of a particle in 
space imposes to add phase factors to the wave function. 
For a massive particle associated with this wave function, 
the two possible phase factors correspond to the energy 
correction e'-^*/^ ~ e""" due to the kinetic energy the 
particle acquires under displacement ; and the displace- 
ment operator of the wave function e'^^/^ ~ ^\mvx/n^ 
thus see from these simple arguments that a full covariant 
phase may include some velocity dependencies. 

Let us quantify this idea. To obtain the complete phase 
transformation, wc will use the covariance properties of 
the Schrodingcr equation under a non-inertial transfor- 
mation defined by x' = x — r {t) and t' = t. Then, the 
differential operators appearing in the Schrodingcr equa- 
tion transform according to d/dt' ~ d/dt + v (t) -V and 
V' = V where v (t) = r is time dependent. Moreover, 
one can show that the Schrodingcr equation for particle of 
mass m and charge q, 



u 



ft2 



*(x,t) , (5) 



is formally equivalent to the one with electromagnetic po- 
tentials V and A, potential U and wave function ^ ex- 
pressed in the displaced (primed) system with the corre- 
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spondence laws 



U' (x',t') 
V (x',t') 
A' (x',t') 



[/(x,t)- 
A(x,t) 



mv (t) 
V- A(x,t) 



(6) 



for displaced potentials and 



(x',i') =exp 



dt 



TOV • X 



(7) 



for the wave function (to being the initial time when 
V (to) = 0). The simplest case of the Galilean trans- 
formation [20| with V = is the obvious limit ^P' = 
exp [i (Et — p • x) /H] of the transformation ©-([Tl (In 
short: Galilean relativity constraints the free particle to 
generate plane waves under displacement). However, for 
treating oscillations we need to consider non-inertial trans- 
formations. Note that the covariance of the Schrodingcr 
equation would be destroyed if r were a function of both 
time and position, because of the appearance of mixing 
derivatives [H]. The assumption (i) above is needed to 
keep the covariance. 

The meaning of the above transformation is clear. 
When a part of a circuit is moving as in Fig|T]&, we have 
two choices. We can either solve the Schrodinger equation 
in the moving frame and to select solutions which satisfy 
the condition of the phase continuity. Alternatively, we 
can transform the equation to the rest frame. In the lat- 
ter case, the phase continuity requirement follows from 
the redefinition ([7]) of the wave function in the displaced 
region. Because the same argument holds with the vec- 
tor potential - we once again have the choice to calcu- 
late the solution when the covariant term A is included in 
the Schrodinger equation and the solution must explicitly 
verify phase continuity after all, or we can simply use the 
gauge covariant phase 7 which guarantees the phase con- 
tinuity everywhere - the complete gauge covariant phase 
will be 



Stt / 27r nie 

$0 / *o e 



r-d\ 



(8) 



where 2me is the Gooper pair mass. Obviously, the veloc- 
ity integral in Eq. (jS]) disappears when the superconduct- 
ing condensate remains at rest along the phase path, and 
thus it connects the displacement velocity to the phase 
difference in the Josephson system and eventually to the 
Josephson current. Note that the electromotive term does 
not require magnetic field. Thus, even a single Josephson 
junction, or a bulk superconductor, may exhibit such an 
electromotive effect (see [T31II1] for more references). 

The inclusion of the relativistic correction in (jS]) sug- 
gests that some subsequent contributions will be found 
for its time derivative. Indeed, the second Josephson rela- 
tion connects the time derivative of the superconducting 
phase with the energy difference across the weak- links, i.e. 
to the electromagnetic work the superconducting charges 
undergo when travelling across the junction. The total 



electromagnetic work is found according to the complete 
Lorentz force. One obtains 



d^p 27r me d 
dt $0 e dt 



vdl 



27r 

$0 
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v X B) -dl 



(9) 



for the phase- voltage relation. For v = 0, we obviously 
recover the (second) Josephson relation ip = 2eV/h. In 
addition to the first term, which generates electromag- 
netic fields due to the motion of accelerating charges, the 
last term in Eq. <Q is another electromotive one, caused 
by the moving surface the magnetic flux threads. Thus, 
even a static magnetic field can generate such a retarded 
effect, because it is due to the kinematic displacement of 
the interfering paths in the presence of an electromagnetic 
field. Indeed, in Ref. [35], it is demonstrated that (for any 
B and S fields) 
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-r / / B dS 



9B 

— V X V X B 

dt 



-v(V-B) 



(10) 

when the infinitesimal surface element dS is time depen- 
dent, with V [t) the velocity field of the contour. From 
Eq.®, Eg. ((To)) reduces to Eq.® using Faraday's law. 
The two expressions (|5|) and © and their explicit deriva- 
tions are the main results of this Letter. We believe that 
they are rather general, because the notation f = J ^ip.dl 
above was a generic phase difference along a superconduct- 
ing path, eventually closed. In particular, it was neither 
ip+ nor iy9_ used for the specific example of the SQUID 
discussed otherwise. 

The equations and (jH]) representing the electromo- 
tive contributions to the constitutive relations provide 
the description of the combined mechanical plus super- 
conducting systems. It is not immediately clear for us 
whether these relations can be easily reproduced starting 
from the phonon concepts, or more generally from micro- 
scopic models of superconductivity (see e.g. [B]). 

Let us now make some remarks about our above deriva- 
tion. From the transformation law (jT]), it is easy to show 



that the electromagnetic sources p 
lm{^*V^} /m transform as (cf. [101 



q\"^\ and '] = qh 



p' (x',t')-p(x,t) 
j'(x',i')-j(x,t)-p(x,t)v(t) 



(11) 



Here, the transformation law for the current is nothing else 
than the London current substitution when a supercon- 
ductor is displaced [7] . We have thus shown that the Lon- 
don current substitution can be demonstrated using the 
covariance of the Schrodingcr equation. We even found 
that the London current is still a correct expression for 
non-inertial displacements. Note that our argument is just 
the quantum version of the London's original one, because 
phase continuity and current conservation are related to 
each other. 

Note also that Eqs.® and ^ are presumably invalid for 
a true velocity field v (x, t). In particular, our assumption 
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(i) was crucial to obtain the London expression for the 
current (recall that j is defined as a gradient interference). 
We also note that the full transformations ([5]), © and 
pT|) are compatible with the usual electromagnetic gauge 
transformations, and that the space-time transformation 
laws for the electromagnetic potentials (and subsequently 
for the electromagnetic fields) correspond to the magnetic- 
Galilean-limit when the electric displacement current does 
not exist, see |2Q1[23| and note [24l. 



Modification of the constitutive relations of the 

SQUID. — We now return to the main discussion and 
include the electromotive contributions (|S]) and © into 
the specific equations describing the SQUID dynamics. 

Usually, the velocity integral term in Eq.® is trans- 
formed to the flux of V X v, then restoring the usual Lon- 
don electromotive effect in term of Josephson current |15) . 
In our case, the integration contour 1 (t) is made of the 
moving bar itself, and thus depends on time as the me- 
chanical parts are moving. Using assumption (i) , we write 
/ v.dl « / (dl) /dt. The integral then reduces to the two 
contact points K12 on Fig|Tj6 and can be approximated by 
£x. Note this is also the rate of surface motion of the cir- 
cuit, or, in our simple geometry, the rate of the length path 
variation. In other words, when a superconducting path is 
moving, the superconducting phase behaves according to 
the static and kinematic properties of the displacement. 

In our SQUID setup, the phase difference thus reads 



k . 

ip^ = —ipe — K-X H X 

Wo 



(12) 



with k K, [meUJo /H) i within our approximation. This 
equation generalizes Eq.®. 

Note that k is the inverse of an effective Compton wave- 
length [TT], related to the inertia of the electrons because 
of their lag behind the lattice mechanical oscillations. Seen 
as a relativistic correction, the k term in (fT2)) can be inter- 
preted as an angular momentum-orbit interaction, in full 
analogy with spin-orbit effects into Josephson physics [53] . 
The main difference between the spin-orbit and the angu- 
lar momentum-orbit interactions is just the chosen frame 
of reference. Indeed, the spin-orbit effect is due to the 
change of frame between the moving electron and the nu- 
cleus at rest in an atom, whereas the angular momentum- 
orbit interaction comes from the change of frame between 
the superconducting charges at rest and the ionic lattice 
which moves with the mechanical vibrations. It is then 
obvious that the mechanical vibrations, like any motion 
of a superconductor, can be seen as the generation of an 
internal angular momentum, which may eventually lead 
to macroscopic magnetic field correction. 

Comparing the two oscillating contributions in (fT2)) 
leads to the ratio B/ujo k m^j^e k, 3ng/C which is very 
small, thus explaining why such a relativistic contribution 
has not been found in previous experiments [31[5]. Never- 
theless, being a velocity term, the Compton contribution 
has to be compared with the Q-factor of the oscillator. 



which may eventually be large enough to make interesting 
relativistic effects observable at the nanoscalc. Also, it is 
interesting to mention that the k term in (jl2p survives in 
the absence of magnetic field. Thus, the observation of 
this electromotive effect in the absence of external mag- 
netic field will be a clear demonstration that a displace- 
ment of charges generates a complete electromagnetic field 
at the quantum level. 

In a SQUID, the phase difference is not affected by 
the contribution ([9]), and the v x B term affects only the 
phase sum iy9+. Eq.® implies that the voltage is defined 
in the following way. 



V 



27r dt 



f+ + K+X X 

Wo 



(13) 



with w K_ « ttB£/^q in our approximation. This 
expression appears natural: because the position of the 
oscillating bar is equivalent to a fiux for the SQUID, its 
velocity i generates voltage. The Compton fc-term plays 
the role of usual Bremsstrahlung in p^ . 

Finally, the above discussions can be summed up in the 
system of equations 



I — sm (^-|_ cos ip- 

RC 
H 

j = sin cos ip^ 



K+i — kx/u)o 



(p+ + K+X 



RC 



Wo 



(14) 



and the constitutive relations 



Lp_ = —ipe ~ K.-X + kxjujQ — f3j 

tp^ = 2ttV/^o ^ i^+i + kx/ujQ 



(15) 



which represent the coupled dynamics of a SQUID with 
an embedded mechanical resonator. The motion of the 
resonator x (t) is governed by Eq.©. 

Self-sustained oscillations. — Now we discuss gen- 
eration of self-sustained oscillations by static current. 

In the overdampcd limit and without self inductance 
{i.e. RC/uJc ^ and /? —5- 0), the macroscopic time aver- 
aged voltage can be approximated as [TH] 




COS^ Pe + K^X 



Wo 



(16) 



in the limit wq/wc <?C 1 (this limit is the correct one for 
actual experiments [HIS], when wq/wc < (lO"'^ — 10^"*)). 
Then, in the oscillating bar equation of motion ([2|), one 
has to realize that the Lorentz force is present only in the 
resistive branch of the SQUID response, and Eq.® can 
be approximated as 



Wo . , 2 (^) 



(17) 



which exhibits self-sustained oscillations when 
gk/RIoujQ > Q~^. Thus, to apply a static current 
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to a SQUID may generate self-sustained oscillations of a 
mcclianical mesoscopic device. This striking effect can 
be observed with longer oscillators than those currently 
used (recall that gk cx or with higher Q-factors. 

Note that the self-sustained oscillations described in p7|) 
originate from the Compton term k. Thus, the presence of 
self-sustained oscillations is independent of the presence 
of the external magnetic flux. 

Conclusion. — In comparison with the initial equa- 
tions of motion (HJ-Q, we have thus established and dis- 
cussed two different kinds of electromotive effects corre- 
sponding to the terms k+ and k in the systems (|14p and 
(fT5|) . These electromotive effects may be relevant for the 
description of the measurement process of a quasi-classic 
elastic bar motion, and/or for the description of the back- 
action effects a SQUID exerts on an embedded mechanical 
oscillator. More precisely, the Compton-like k term seems 
to be important for backaction effect, because it generates 
magnetic flux ; whereas the kinetic interference term 
transforms existing magnetic field into voltage, and might 
be important for detection purpose when the voltage is 
monitored. 

We demonstrate that under certain conditions, the elec- 
tromotive effects can lead to the development of self- 
sustained oscillations in the suspended SQUID geometry. 
This is, however, not the only explanation: Capacitive 
effects may be responsible as well for self-sustained os- 
cillations, without requiring an explicit velocity depen- 
dent phase-flux relation {i.e. in the limit kx/uj^ — and 
K+x — in (|15|) but with a finite RC /lUc term for which 
i = {V) /Rio + Cd{V) /lodt with the Eq.^ as a trial 
expression for {V) and with k ^ 0). One possible way 
to distinguish between these two electromagnetic contri- 
butions is to consider vanishing magnetic field, see [HI for 
more details. 

In order to find the relativistic contributions, we as- 
sumed that it is formally possible to separate the two sub- 
systems, and use some relativistic (i.e. covariance) argu- 
ments of their quantum properties. In addition, we believe 
that such electromotive contributions are difficult to ob- 
tain using phonon models, which opens the question of 
appropriate arguments to describe elasticity and mechan- 
ics at a mesoscopic length scale, in particular the inclusion 
of mechanical systems into electromagnetic ones. Testing 
our conclusions beyond these assumptions goes beyond the 
scope of this paper. 
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